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Abstract. In quantum/wave systems with chaotic classical analogs, wavefunctions 
evolve in highly complex, yet deterministic ways. A slight perturbation of the system, 
though, will cause the evolution to diverge from its original behavior increasingly with 
time. This divergence can be measured by the fidelity, which is defined as the squared 
overlap of the two time evolved states. For chaotic systems, two main decay regimes of 
either Gaussian or exponential behavior have been identified depending on the strength 
of the perturbation. For perturbation strengths intermediate between the two regimes, 
the fidelity displays both forms of decay. By applying a complementary combination of 
random matrix and semiclassical theory, a uniform approximation can be derived that 
covers the full range of perturbation strengths. The time dependence is entirely fixed 
by the density of states and the so-called transition parameter, which can be related 
to the phase space volume of the system and the classical action diffusion constant, 
respectively. The accuracy of the approximations are illustrated with the standard 
map. 
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1. Introduction 

In our contribution to this special issue on random matrix theory (RMT), we address 
the subject of wave field evolution governed by linear wave equations. In particular, 
our interest is in the sensitivity of evolution to perturbations. The treatment of this 
problem illustrates two well-known, important lessons arising in RMT studies. First, 
RMT and semiclassical theories, when used in tandem, give a much stronger theoretical 
approach to a problem than either was capable of on its own [Q. Second, even though 
RMT is inherently unphysical in certain ways, nevertheless, careful application leads to 
physically meaningful predictions |^. For example, using RMT to mock up the evolution 
of a wave packet generates a dispersal into the entire Hilbert space simultaneously. 
Thus, it lacks a dynamical continuity property that a semiclassical theory would easily 
reveal as absurd. It turns out that for a wave mechanical system with a chaotic 
underlying ray dynamics, there is a logarithmic time scale dependence on the wave 
vector p, ^, |], and RMT does not "know" about this; i.e. its results cannot be 
trusted for wave propagation up to this time scale. Yet, beyond the so-called logtime 
(where the semiclassical dynamics increases in complexity exponentially 0), RMT can 
be a powerful tool for understanding the statistical nature of the dynamics. In what 
follows, we take advantage of the RMT-semiclassical complementarity and avoid making 
inferences about pre-logtime dynamics. 

The study of the sensitivity of wave field evolution to perturbation has recently 



attracted a great deal of interest |0, ^, ^ [T0|, |ri|, |T2|, for several fundamental reasons, 
i.e. decoherence, reversibility, tomography... In the case in which the basic system is 
weakly coupled to an environment, the sensitivity relates to the manner in which the 



system decoheres [|Tj]. The recent focus on quantum computation has greatly increased 



the attention given to this issue |T^. Another motivation has been the reversibility 



of wave field evolutions such as the spin polarization echoes in nuclear magnetic 
resonance |T^. For a broad class of problems, propagating forward in time with two 
slightly different systems is mathematically equivalent to propagating forward with one 
system and reversing the dynamics with another for the same propagation time. Another 
application is to problems of waves propagating through random media (WPRM), 
the medium is likely partially characterized and often slowly evolving; twinkling of 
starlight passing through the atmosphere is an example. The interest is often in 
performing tomography or at least knowing what information can be deduced about 
the medium |T^. Similar statements can be made regarding the study of disordered 



systems [0. Lastly, we mention that problems involving simple, chaotic systems or 
even strongly interacting, many-body systems usually require taking measurements as a 
function of some system parameter. For these systems, single measurements may almost 
seem like random number generation, whereas it turns out that important information 



can be extracted from the parameter variations of the measurements |]I9[. The trick is 
to know what information therein exists. Our focus here will remain on the wave field 
evolution behavior without getting into specific physical applications. 
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Introduced nearly twenty years ago, a natural measure of the sensitivity of wave 
field evolution is the fidelity, which is defined as the squared overlap of some initial state 
propagated forward in time with two slightly different systems |^|. For normalized 
wavefunctions, it begins at unity and typically decays to some small value. In the spin 
polarization experiments previously cited, it was observed that a change in the functional 
form of the fidelity's decay occurred depending upon the dipolar interactions The 
decay exhibited either Gaussian or exponential behaviors and in some instances both. 
This mirrors known behavior found in simple, quantized chaotic systems P, note 
that there also exists some work on integrable, near-integrable, and mixed phase space 
systems |2^, ^ which exhibit other behaviors. For extremely weak perturbation 
strengths, quantum perturbation theory along with some simplifying assumptions 
quickly lead to the prediction of a Gaussian decay [Q. For larger perturbation 



strengths, this approach fails, and a Fermi Golden Rule regime emerges which implies 
an exponential decay behavior p5[. Semiclassical theory gives this decay rate in terms 
of a classical action diffusion constant . If the perturbation strength is even larger, at 
some point the fidelity decay rate saturates because it is impossible to decay faster than 
the logtime this is known as the Lyapunov regime because the rate is independent 
of the perturbation strength, and often determined by the Lyapunov exponent. 

Regardless of the strength of the perturbation or the nature of the dynamics, for 
very short times the fidelity is quadratic. This regime which can be derived by using 
time dependent perturbation theory may be extremely short in time and difficult to 
observe. Likewise, in the long time limit, the fidelity will saturate to the inverse of the 
effective number of states in the system. Both of these effects decrease in importance in 
the asymptotic limit, and are disregarded in what follows; the Lyapunov regime is also 
disregarded since that involves pre-logtime dynamics. 

In this article we concentrate on developing a uniform approximation that 
encompasses both the Gaussian and exponential decays as well as the intermediate 
regime for simple chaotic systems; and we restrict our attention to the Schrodinger wave 
equation. The fidelity undergoes a crossover from exponential to Gaussian behavior as 
a function of time, and the value of the perturbation strength determines whether the 
crossover is visible. In fact, it turns out that the only information residing in the decay 
of the fidelity is the density of states and the transition parameter of RMT |^|; the 
spreading width could also substitute for the transition parameter. Alternatively, using 
semiclassical arguments, one could also say that the only information contained is the 
phase space volume and the classical action diffusion constant 0. 

The paper is organized as follows. In the next section we review the Gaussian and 
exponential Fermi Golden Rule regimes found in simple quantized chaotic systems. In 
Sect. 3, a uniform approximation is derived using RMT and a degenerate perturbation 
theory. The theory is applied to the quantized standard map in Sect. 4. It is found to 
give an excellent approximation of the quantum results. We then conclude with some 
discussion and closing remarks. 
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2. The Fidelity 

Consider an evolving quantum wave function \ax{t)) where A defines the parameters of 
a Hamiltonian, H{X). The overlap of the same initial state, \a), propagated via two 
different Hamiltonians specified by Ai and A2 is 

Ax{e-t) = (ttAiWIttA^W) = {a\Ul{t)UxMc^) (1) 

where Ux{t) is the unitary evolution operator, e = Ai — A2, and A = (A2 + Ai)/2. In 
the limit that the system is strongly chaotic, the statistical properties of the overlap 
are independent of A, so we omit the A subscript from here on. The fidelity is just the 
absolute square overlap 

C{e;t) = \A{e;t)\' (2) 

It is sufficient to work with the overlap itself in developing the theory since the operation 
of squaring only introduces cross- correlation effects higher order in e that can be ignored. 

By inserting complete sets of states, the overlap can be rewritten in terms of the 
eigenenergies and eigenstates of the perturbed and unperturbed systems as 

= im("(0)l'^Ai)(nAjmA2)(mA2|a(0)) 

n m 

X exp{-i[E„(A2) - E„(Ai)]t//i} (3) 

The Bohigas-Giannoni-Schmit conjecture |2^ and later works suggest that for a chaotic 
system, it is appropriate to apply RMT arguments to understand its statistical behavior. 
Here, as mentioned already, one has to add the caveat - as long as e is weak enough 
that C(e; t) decays on a much longer time scale than the logtime. The first property 
we would like to invoke is the invariance of the ensemble over changes in bases. Thus, 
any rotation to a different initial state must generate the same decay behavior for the 
fidelity. It suffices to make an average over a complete set of initial wave functions |a(0)) 
to obtain 

= ^ E E I (^A, |mA,) r exp{-i[E„(A2) - E^{X,)]t/h} (4) 

n m 

where is the number of states in the complete set. RMT has been shown to be a 
strongly ergodic theory 0, which in this case implies that the fluctuations of the 
quantity 6A = A{e; t) — A{e; t) vanish with increasing matrix dimensionality as A^~^. 
The averaged overlap may be decomposed into a diagonal and off-diagonal part 

= ^EK^AjnA,)pexp[-iAE„t/;i] 



N 
1 

N 



+ 4E E \{nxAm,,)\'exp{-i[E^{X,)-E^^{\,)]t/h} (5) 

where AEn = En{\2) — E„{Xi) is the change in an eigenenergy due to the perturbation. 
For e = 0, the off-diagonal term exactly vanishes, the diagonal term equals unity, and 
the fidelity does not decay as must be case. 
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The quantum perturbative regime 

The diagonal term is defined by adiabatically following an eigenlevel as the system's 
parameters are continuously varied. It may be further evaluated by ensemble averaging 



1 



n 



-iAEnt 



(6) 



h 

It is most significant for extremely small perturbation strengths. From quantum 
perturbation theory, the lowest order shift of an eigenvalue depends on the diagonal 
element of the perturbation whereas the rotation of the eigenstates depends on off- 
diagonal elements and energy differences. Within RMT the perturbation matrix 
elements are uncorrelated and thus to lowest order in e, the ensemble averaging of the 
amplitudes and phases can be separated. The average over the phases is over a Gaussian 
probability density for the diagonal matrix elements. Also, the ensemble average of the 
amplitude is independent of n. Thus, 

(I^-RT)) ^ (|KJn,,)p)exp(-eV,¥/2;.2) (7) 

where cr^ is the variance of the diagonal matrix elements of the perturbation V defined 
by writing the Hamiltonian locally as H{X) = Hq + eV. In the limit of a differential e, 
is equal to the variance of the level slopes by the Hellmann-Feynman theorem . If 



the perturbations are sufficiently small such that the eigenfunctions do not significantly 
vary, i.e. {nx^\mx2) ~ Snm, then the off-diagonal terms will be negligible and may be 
ignored. Hence, the fidelity has a Gaussian decay in time in the weak perturbation 
regime. 

The diagonal matrix element (level slope) variance is not a free parameter, and can 



be given an elegant interpretation through a semiclassical approach I^O], |3T|, It is 
given by 



2 



0-,, 



2gK{E) 



(8) 



where 2g/(3 is the number of classical orbits with identical action and d is the total 
mean density of states. The index (3 = 1 for time-reversal-invariant systems and j3 = 2 
for time-reversal-breaking systems. K{E) is the classical action diffusion constant on 
the energy surface E given by 



K{E)= {V{p{0),q{0);X)V{p{t),qit);X))^Jt (9) 



In this expression the averaging is defined over the primitive periodic orbits of long 
period. The physical picture is that the action difference of two long orbits continuously 
deformable into each other as a function of e is the result of a diffusive process. As the 
orbit explores the phase space "randomly", sometimes its action increases relative to 
the other orbit, sometimes it decreases. The variance of that diffusion is proportional 
to the variance of the level slopes. 
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Semiclassical analysis 



The above analysis describes a very restricted range of perturbation e centered at zero. A 
semiclassical approach without the quantum perturbation limitations can be developed 
that is valid over a much broader range. The semiclassical construction of an evolving 
wave function begins with the propagator 



(q|f>|q') 



1 \ 



d/2 



2mh 



) S 



det 



'9W,(q,q';t)' 



1/2 



9q9q' 
X exp ( ipy, (q, q'; t)/h 



(10) 



The phase is specified by the time integral of the Lagrangian Wj{q,ci';t) and and an 
index based on the properties of the conjugate points (like focal points), Uj. 
The overlap decay in terms of the propagator is 

A{e;t) = I dqdq'dq''(«|q)(q|f>ljq')(q'|f>Ajq")(q"l«) 



For the reasons given in the discussion surrounding equation (^ , there can be no overlap 
decay dependence on the specific initial state for chaotic systems as long as the decay is 
on a time scale which is much longer than the logtime. Since we have explicitly excluded 
the Lyapunov regime from consideration in this paper, we are free to choose any form 
for the initial state. For convenience, we take Gaussian wave packets 

\2 



(qla) = (tto"^ 



-d/4 



exp 



iq- qo 



2(t2 



+ ^(q-qo 



(12) 



The initial coordinate of the actions are expanded about the center of the wave packet 

w^i(q, q; t) = w^i(q, q^; ^) - Pa ■ (q' - qa) + 

.T 9iy,(q,q';t) 



2^^ 



qo 



iq - qo 



9q'9q' ""^ ^^^^ 

Changes in the amplitudes are ignored, since the most important contribution to the 
overlap decay is each orbit's action change due to its division by h,. Applying stationary 
phase on the q and q" integrals, we obtain 



A{e,t)^ jY^Cjexp ^AVrj(q, q«; t) 



dq 



(14) 



where Cj is a magnitude which disappears from the expressions ahead, but could be 
deduced using the equations in the appendix of [0. The AWj{q, q^; t) are the action 
differences of two rays that begin at qo, and end at position q in a time t and are 
continuously deformable into each other; bifurcations are neglected. First order classical 
perturbation theory gives 



AW^i(q,q«;^) =e — —g^ dt 



(15) 



The stability matrix elements and the Lagrangian are evaluated along the orbits of 
H{X). In general, stationary phase integration cannot be performed on the last integral 
because the action differences are less than h. 
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Again a statistical argument can be employed. The weighted phases will be nearly 
Gaussian distributed for the reason mentioned earlier that long orbit action changes are 
part of a diffusion process 0. The last integral and summation in equation (14) can be 
replaced by a Gaussian integral. Thus, we obtain 

A{e;t)^exp{-a^w/2h^) (16) 

where a^y is the variance of the action differences and is given by pO| , pT 



cr^ = 2e^K{E)t = {2'Kdal)e'^n(3t/2g (17) 

with the same action diffusion constant as before. In contrast to the quantum 
perturbative argument giving a quadratic t-dependence in the exponential, here the 
argument is linear. The semiclassical argument does not contain the action correlations 
necessary for the theory to contain the quantization of the spectrum. It therefore 
misses the Gaussian contribution from the diagonal terms, and essentially evaluates 
the off-diagonal contribution in the regime of larger perturbations as though there is a 
continuum. The last form given in the equation exhibits the Fermi Golden Rule form 
expected for such circumstances. Note that the symmetry factors do not enter K{E) 
(the symmetry dependences cancel in the final form) in contrast to the diagonal term. 

3. A Uniform Approximation 

A straightforward argument determines the value of e which separates the exponential 
and Gaussian regimes. We emphasize that this is not the cross-over to exponential 
dependence beginning from the quadratic time dependence derived using time- 
dependent perturbation theory, but rather the relative strength of contributions from 
diagonal versus off-diagonal contributions in equation (||). Note that either dependence 
dominates over the time range in which the other's argument of the exponential takes 
on the lesser value. A cross-over between the two regimes occurs when their arguments 
are equal. Simple algebra gives the this time scale t* = hdf3/2g, which is the Heisenberg 
time, th, to within the symmetry factor. Exponential decay dominates if the decay is 
completed by th, and Gaussian decay dominates if little decay has occurred by th- In 
between it turns out that the decay has components of both behaviors. In terms of the 
parameters, the question is whether 



hg 



K{E)t* Txd(3K{E) ' ' 

It turns out using the results derived next that near the equality, the magnitude of the 
off-diagonal term is comparable to the diagonal term. 

Returning to equation (0), it is an oversimplification to take the amplitude equal 
to unity. Instead, we define the fraction 

/ = (K^A,K)p) (19) 

and ahead give a theory for its value in the neighborhood of unity (small e) using a 
degenerate perturbation theory and RMT. It is worth noting the broader significance 



Uniform Approximation for the Fidelity 



8 



of the quantity /. It and some closely related measures have previously been 
studied by several authors. In ||3^, the authors give the leading analytic correction 
from unity for the case of broken time reversal symmetry, and infer the asymptotic 



behavior for large e, and in [^S], |36|, a rough guess is given for the full range of 
behavior consistent with numerical simulations with random matrices. These measures' 
interpretations as parametric correlators of eigenstates, eigenstate components, or 
widths were emphasized. Motivations for their introduction come from a broad variety 
of possibilities of parametric dependence on controllable external or uncontrollable 
variables such as electromagnetic fields, thermodynamics quantities, or shape/geometric 
properties. Physical realizations from nuclear fission to quantum dots to atomic and 
molecular spectroscopy are possible. It is interesting to note that the same quantity is 
entering into the theory of the fidelity, but only for transition values of e between the 
two domains giving exponential and Gaussian decay. 

The argument giving a Gaussian behavior for the phase average remains valid to a 
regime where / is quite small and the diagonal term negligible. It suffices to evaluate / 
to complete a more broadly applicable expression for the diagonal term. 

The off-diagonal term is more subtle. By normalization, the prefactor must be 
equal to 1 — / , but the argument of the exponential needs to be determined. Recall that 
semiclassical dynamics is valid far beyond the logtime, but breaks down on algebraic 
scales before th 0, |3^. In particular, this implies that in the limit of a very weak 
perturbation, the exponential scale as previously derived places the exponential behavior 
beyond th- The semiclassical expression turns out to require modification. Since it is 
valid on time scales shorter than th, we can expand the exponential to the linear time 
term. This implies that in this regime the argument of the exponential multiplied by 
1 — / must be equal to the argument found in the previous subsection. Therefore, to a 
good approximation, the overlap can be written 

— — — ^ „ -e'^a^t^ ,^ -e^Ti(3aldt 
A(., t) « /exp + (1 - /) exp ^^-^ 

where the first form emphasizes the quantum information (matrix element variance or 
level slope variance), and the second form emphasizes the classical information (classical 
action diffusion constant). The transition between the two regimes, Gaussian and 
exponential, as the strength of the perturbation is increased, is given by the fraction / 
remaining in the diagonal overlap factor. The exponential decay rate has been modified 
by the introduction of 1 — / in the denominator which is necessary to offset the reduced 
normalization. In the limit of large perturbations where / — > we recover the previous 
results for the exponential Fermi Golden Rule decay, and in the small perturbation limit, 
/ — s> 1, the Gaussian decay. 
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Degenerate perturbation theory within RMT 

To derive the leading corrections to / for both time-reversal-invariant and noninvariant 
cases, the most direct approach would be to use the second order perturbation theory- 
expression 

where Vmn = (^aJ^I^Ai)- Within RMT, is a Gaussian random matrix with the 
correct symmetries, Gaussian orthogonal ensemble (GOE), unitary ensemble (GUE), 
etc... We only consider the time- reversal-invariant and time- reversal- noninvariant cases, 
GOE and GUE respectively. 

Define an energy variable s rescaled by the mean local spacing and a unit variance 
variable x for the perturbation matrix elements. Then the above sum can be reexpressed 
as the following double integral 



oo 



/ ~ 1 _ 2A / / ^R^(s')pix; s')ds'dx (22) 



s'^ 



where A = t^v'^d^ is the transition parameter and v"^ is the local variance of the 
off-diagonal matrix elements of V p6i. d^ and are defined for the irreducible 



representations of the system. If there are symmetries, then these quantity must be 
adjusted accordingly by dividing out these symmetries, eg. dr = d/g and = (3al/2. 
The function R2{s) is a sequence of (5-functions which determine the spacings, and is 
known as the two-point density correlation function. Similarly, the function p(x; s) is 
a sequence of 5-functions whose intensities determine the matrix elements. Crudely 
speaking, taking the ensemble average of the above expression is reduced to inserting 
the R2{s) function from the appropriate Gaussian ensemble, and taking p{x; s) to be a 
zero-centered, unit variance, Gaussian probability density independent of s. 

It was emphasized by French and co-workers that the transition parameter A 
was quite general Although it naturally emerged from perturbation theories, its 
relevance extended to entire symmetry breaking transitions and it is the only scale that 



can enter those problems; a similar statement applies to parametric statistics as well |38 

It is straightforward to understand that for s near zero, the integral diverges for 
the GOE, and that this approach must fail to give proper results for /. The two- 
point correlation function is approximately linear, R2{s) ~ n'^s/G, so the above integral 
will diverge logarithmically. However, we can apply a degenerate perturbation theory 
approach developed in |^ ^ ; the degenerate theory also gives an extra valid correction 
for the GUE. By diagonalizing a two by two matrix, the difference in eigenenergies is 
given by |0 

1 



^ m 

X sign(E„(Ai) - E™(Ai)) (23) 
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Using this result, equation (E2|) can be reexpressed as 



Away from very small s', the degenerate perturbation result is unnecessary, but still 
gives the exact leading order results as a function of A. This integral is rather difficult 
to evaluate exactly, but the correct leading order results can be extracted by splitting 
the integral into two regions [0, s] and [s, oo\. First, the eigenstate overlap is written 



/ 



/°c> 2A 1 



(25) 



^2(r) = r"' . L (26) 



where the power series for R2{s) is substituted in the first term, and the exact function 
is left in the second term. For our purposes, it is sufficient to replace R2{s) in the first 
region by vr^s/G, which gives the GOE result, or by tt^s^/S, which gives the GUE result. 
The exact two-point correlation functions are 

1 - [s(r)]2 + Js{r)Ds{r) GOE 
1 - [s(r)]2 GUE 

where s(r) = sin(7rr)/7rr, Ds{r) = ds(r)/dr and Js{r) = Jq s{t)dt — 1/2. Next, 
the [0, s] integrals are evaluated exactly, and their asymptotic expansions kept up to 
(not including) 0{A'^). The other integral is also exactly integrated, and its small-s 
expansion is made. When properly formulated and the two evaluations summed, all the 
divergent terms less than 0{A'^) vanish; also the s-dependent terms coming from the 
[0, s] integration vanish. The s-independent terms give the desired result. The remaining 
s-dependent terms arising from the [s, oo] integral vanish as higher order terms in the 
power series of R2{s) are included in the first integral. In this way, the final result 
does not depend on the value chosen for s. Carrying out this technique up to but not 
including ©(A^), we obtain 



1 - 4^ ( 1 - 7 - In 



GOE 



f ^ { , 2n^A I 327rA /?;— T ^ttt^ (2^) 



, . 3 +^V27rA GUE 

where 7 is the Euler constant. The GOE results and the correction term (A^^^) for the 
GUE are new, and the leading GUE term is consistent with reference [Q, which also 
gives the asymptotic results 

f^l^ (28) 
GUE ^ ^ 

The overlap ensemble average is expected to be universal and should apply to 
all systems chaotic enough to expect RMT behavior. Numerical studies and heuristic 
arguments suggest that the overlap between perturbed and unperturbed eigenstates is 
roughly a Lorentzian-like function of the perturbation strength [^, We look for a 
dependence of the form 

^ = Tt1(a) 
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Note that in the previously discussed forms, the expression for / was raised to the power 
f3, h{A) ^ apA, and approximate values of the constants were given by a numerical 
fit. Using perturbation theory Wilkinson and Walker [|T| gave an expression for the 
individual off-diagonal matrix elements which has a compatible form. Nevertheless, the 
O(AlnA) and 0(A^/^) corrections plus the asymptotic expressions indicate difficulties 
with their interpolation since it gets both limits of small and large A incorrectly. 
Equations (|27| - |29|) are, of course, insufficient to determine h{A). However, we suggest 
the following as perhaps the simplest forms for somewhat improved interpolations that 
are at least analytically consistent with the limiting regimes of large and small A: 

TT^A (t _ ^ _ ( l + boTT^A^ \ 1 TT^A \ 

6 V ^ Vl+fciTrSA+boTr^A^y 4+exp(5+7)7r2A ^ VjWI^ 

K^)^{ „ , , (30) 



27r2A ( l+fcoTT^A 



GUE 



3 

The use of a quotient of two polynomials inserts a general function with some fitting 
parameters that we determined with GOE and GUE Monte Carlo simulations; see the 
appendix for the Monte Carlo details. Our best fits are: 6o = 0.932, hi = 0.618, 
and Co = 1.89. Figure |l] shows the differences between the Monte Carlo results and 
the new interpolations along with the differences between the Monte Carlo results and 
those from references Even though the new interpolations have the correct 



limiting forms, they still appear to have some small, detectable deviations from the exact 
RMT predictions. Perhaps, the main conclusion is that the previous interpolations do 
surprisingly well given that they fail to get any of the limiting cases correctly. 

Again, semiclassical theory fixes the value of A using the arguments presented in 
the previous section. Thus, 

e^^_eW(E) 
^" 2g^ ~ Tigh ^^^^ 
where the last form follows from equation (||). In some simple systems such as the 
standard map used in the next section, K{E) can be calculated analytically. In fact, 
the criterium, equation (plSj) , for observing either Gaussian or exponential behavior can 
be written in terms of the transition parameter 

tt'A > i (32) 

It is remarkable that the same transition parameter that determines the effects of weak 
symmetry breaking in RMT, gives precisely the exponential-to-Gaussian cross-over for 
the fidelity. 

Finally, the uniform approximation for the fidelity is the square of equation (pO|) 
with the results given in equations ( P^ - PT|) . It could have turned out that there were 
three parameters in the functional form of the decay, the relative amount of Gaussian 
to exponential component (/), and the scales in the arguments of the Gaussian and 
exponential functions. Interestingly, other than the density of states (or phase space 
volume), only one quantity determines all three potential parameters, that being the 
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classical action diffusion constant. There is no other information contained in the decay 
of the fidelity for a strongly chaotic system. 

4. Fidelity in the quantized standard map 

We use the standard map, which is a paradigm for chaotic systems, to demonstrate our 
results for time-reversal-invariant systems. For a kicking strength A greater than 6 or 
so the standard map is fully chaotic. The classical map is defined by 

Pi+i = Pi - (A/27r) sin(27rgi) mod(l) 

qi+i=qi+Pi+i mod(l) (33) 

The quantized propagator with discrete levels is given by 

{n'\U\n) = r — exp[i7r(n — n')^/A^] 
yiiV 

/ y^jY \ 
X exp I i— cos[27r(n' + a) jN] J (34) 

where n^n! = 0, . . . , N — 1 and a is a phase term which we set equal to zero. The 
effective Planck constant is h = 1/N and the average density of states is d = N/2tt. 

From the considerations of Sect. 2, examples of nearly pure Gaussian or exponential 
decay can be found by selecting an appropriate value of e using equation (p^). In figure 
we see that for small h (large A^), even without initial state averaging, the theoretical 
curves for the Gaussian and exponential regimes respectively match the quantum results. 
The theory curves use the analytic form of the action diffusion constant, which for the 
standard map is [pO| 

K{E)^[l + 2J2{\)]/A{27iY (35) 

Also, it turns out that the expressions for the perturbation matrix element variance 
or level slope variance is different for quantized maps than for continuous dynamical 
systems since the Floquet eigenangles are scaled differently than energy eigenvalues. 
Here, the variance of the eigenangle slopes is 

al = al/h^ = g{ATifNK{E) where g = 2 (3 = 1 (36) 

In figure ^ we show an example of the behavior of / as a function of the transition 
parameter A. Ergodicity allows the ensemble averaging to be replaced by an average 
over eigenstates. For quantum maps, the transition parameter is defined in terms of the 
variance of the eigenphases as opposed to the eigenenergies as in continuous systems. 
Therefore, for the quantized standard map 

A=—^ = — I + 2J2A 37 

2^2 g 8(27r)4^ ^ ^ ' 

Using the GOE interpolation expression for / in equation (^) gives excellent agreement 
with the quantum results. 
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Finally, we show an example in figure ^ where e has been chosen so that A{e; t) 
exhibits both exponential and Gaussian behaviors in different time regimes. The 
theoretical expression given in equation (^) with equation ([35| ) and the interpolation 
approximation for / of the previous section compares well over the full time dependence. 
It is possible to isolate both the exponential and Gaussian components of A{e;t), and 
these are shown in the lower panels. The intercepts of the curves at t = show the 
accuracy of the numerical fit to /, and the decay scales can be seen to be given properly 
as well. 

5. Conclusions 

The fidelity is extremely useful for studies of decoherence, reversibility, and tomography 
in several contexts, and has received a great deal of attention recently. It is to be 
expected that RMT theory would be helpful in understanding its behavior for strongly 
chaotic systems. Indeed, we applied RMT to develop a uniform approximation that 
covers both previously known regimes of Gaussian and exponential decay, and the 
theoretical predictions describe the behavior of the quantized standard map extremely 
well. Although not noted earlier in the text, the range of A over which the bare 
degenerate perturbation theory is accurate, is disappointingly narrow. However, it 
helped to provide a schema for an improved interpolation formula. We further included 
the semiclassical analysis required to show that only one parameter (other than the 
density of states/phase space volume) determined the full behavior, the classical action 
diffusion constant. In quantum terms, the sole information is carried by the transition 
parameter or alternatively by the perturbation's spreading width. As a function of 
this parameter, the uniform approximation gives a continuous family of universal decay 
curves for all strongly chaotic systems. 

Our numerical tests of the time-reversal-invariant uniform approximation with 
the quantum kicked rotor relied on large dimensionality examples, i.e. = 1000, to 
show that although we averaged over initial states to generate the theory, the ergodic 
properties of RMT imply that the results apply individually to nearly any initial state of 
a strongly chaotic system. For smaller values of A^, the fluctuations due to not making 
an initial state averaging increase and are readily visible. There averaging would be 
required to have the same high quality comparison between the theoretical prediction 
and the quantum results. We would expect a similar quality of results to apply for 
time-reversal-noninvariant chaotic systems. 

The new interpolation formulae given in equation ( pO| ) for the overlap of a perturbed 
and an unperturbed eigenstate are improvements over the simple Lorentzian forms since 
they contain the correct limiting forms, and are more faithful to the Monte Carlo 
simulations. Nevertheless, the previous forms work surprisingly well; even where the 
slope at small perturbation goes to infinity in the GOE results, the previous form never 
goes too far astray. The overlap / is important in other contexts as well. For example, 
in mesoscopic physics / is related to the conductance of a quantum dot undergoing a 
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parametric change. For reasons such as this, it would be interesting to know the full 
analytic transition with perturbation strength between unity and zero. 

For integrable, near-integrable, and mixed phase space systems, the behavior of the 
fidelity is not universal and depends upon the initial state chosen. Similarly, chaotic 
systems in the Lyapunov regime (not studied in this paper) could also exhibit initial 
state dependent phenomena. Such systems are also of interest and contain a variety 



of different regimes; see the work in P2| , P3| . We are currently pursuing several open 
questions in these dynamical systems. 
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Appendix A. Monte Carlo simulations of / for the GOE and GUE 

We use the Hamiltonian 

H = sm{€)Ho + cos{e)Hi (A.l) 

in the Monte Carlo simulations. Hq and Hi are random (300 x 300) matrices in either the 
GOE or GUE symmetry. New random matrices are recalculated for each e to generate 
independent samples. An ensemble of 300 such matrices are used where the average 
density of states, d, and variance of the level velocities, cr|, are averaged over the entire 
ensemble. Only the middle third of the spectrum is used. The eigenfunctions of Hq are 
overlapped with those H and then the absolute value square is taken to obtain /. A is 
given by 

r eV|//2 GOE 
A = _ (A.2) 

[ e^ald" GUE 
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Figure 1. Overlap of eigenstates, /. The two upper and two lower plots are for 
GOE and GUE statistics, respectively. The solid curves in the upper and lower middle 
plots are results from RMT calculations for an ensemble of 300 random (300 x 300) 
matrices. The dashed curves are the uniform interpolation, while the dotted curves 
are the lorentzian ansatz of reference | ]35[ |. The upper middle and lower plots show 
the difference between RMT and the uniform interpolation (solid) and the difference 
between RMT and the lorentzian ansatz (dashed). 



Figure 2. Example of Gaussian and exponential decays. The quantum standard 
map curves are solid and the theoretical curves are dashed. The upper plot is for 
e = 2 X 10^"^ and the lower plot is for e = 10^^. The other parameters are N — 1000 
and A = 18. 



Figure 3. Overlap of eigenstates, /, for the quantum standard map. The solid curve 
is for the quantum standard map averaged over eigenstates for N = 1000 and A = 18. 
The dashed curve is the interpolation formula using RMT and degenerate pertubation 
theory. 



Figure 4. Crossover decay. The upper panel demonstrates the agreement between 
the quantum standard map (solid curve) and the theoretical results (dashed curve). 
The middle and lower panel shows each piece of the fidelity (diagonal and off-diagonal) 
separately and the corresponding theoretical results (Gaussian and exponential). The 
parameters are N = 1000, A = 18 and e = 5 x lO^''. 
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